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ABSTRACT 

Bright and nearby (22 pc) solar-type dwarf k For (HIP 11072) is a triple system. The close pair of 
M-type dwarfs Ba,Bb with a tentative period of 3.7 days moves around the main component A on a 
26-year orbit. The mass of the "dark companion" Ba+Bb is comparable to the mass of A, causing 
large motion of the photo-center. The combined spectro-interferometric orbit of AB is derived and 
the relative photometry of the components A and B is given. A weak signature of Ba and Bb is 
detected in the high-resolution spectra by cross-correlation and by variable emission in the Bahlmer 
hydrogen lines. The activity of the M-dwarfs, manifested by a previously detected radio-flare, is likely 
maintained by synchronisation with their tight orbit. We discuss the frequency of similar hidden triple 
systems, methods of their detection, and the implications for multiple-star statistics. 
Keywords: stars: binaries - stars: individual: HD 14802 
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1. INTRODUCTION 



In January 1993 iGuedel et al. I (Q~995) detected mi- 
crowave radiation from the nearby solar-type dwarf 
k For. Its source was located 0'.'23 south of the star. At 
that time it was already known that k For is a spectro- 
scopic binary. So, the authors suggested that the radio 
emission comes from the secondary companion, presum- 
ably a low-mass flaring dwarf. As we show below, the 
secondary is in fact a tight pair of M dwarfs. Despite the 
old age of this triple system (4 to 6 Gyr) , fast axial rota- 
tion of the M-dwarfs (hence high activity) is maintained 
by synchronization with the orbit. 

Kappa Fornaci (HIP 11072, HD 14802, HR 895, GJ 97, 
a 20 oo =2:22:32.54, <5 20 oo = -23:48:58.8) is located at a 
distance of 22 pc from the Sun ac cording to the origi- 
nal Hipparcos catalog (|ESA 1119971 ) and its new reduc- 
tion (jvan Leeuwen 1 120071 ). The Hipparcos satellite de- 
tected a large acceleration of 19.4 mas yr -2 presum- 
ably caused by the invisi ble (astrometric) companion 
(jMakarov fc Kaplan I [20051 ). iGontcharov et al. I (|2000f ) 
collected historical astrometric data spanning half a cen- 
tury and sugg ested that the orbital period is around 
26.5 yr. Later iGoncharov fc Kivaeva I (|2002l ) published 
the elements of the astrometric orbit with this period 
and with a relatively large semi-major axis of a — 0'.'26, 
noting that the "dark companion" is massive. 

Radial velocity (RV) of k For was monitored both 
at Lick and at La Silla in search of planetary compan- 
ions. None was found so far, but the R V trend caused 
by th e stellar companion was obvious. iNidever et al. I 
(2002) noted the trend, while lEndl et al. I (|2002l ) pub- 
lishcd a preliminary spec t rosco pic orbit with a 21yr pe- 
riod. lAbt fc Willmarthl ()2006[ ) added their own obser- 
vations and revised the orbit slightly, although the full 
orbital cycle was not yet covered. Neither of these pub- 
lications mentions the astrometric results. 
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1 Based on observations obtained with CHIRON spectrometer 
at the 1.5-m CTIO telescope operated by SMARTS (NOAO pro- 
gram 2012B-0075), at the SOAR telescope, and at the Gemini 
Observatory (program GS-2012B-Q-71, PI M. Hartung) 



The astro metric and spectroscopic companion was first 
resolved by lLafreniere et al. 1 (120071) in 2005 a nd, inde- 
pendently, bv iTokovinin fc Cantarutti I (|2008[) in 2007. 
For this reason it received two confusing "discoverer" 
codes LAF 2 7 and TOK 40 in t he Washington Double 
Star Catalog (jMason et al. II2001| )FI By using all resolved 
measures and fixing the orbital p eriod and eccentr i city to 
those of the Gontcharov's orbit, lHartkopf et al. I ()2012l ) 
computed the first visual orbit of grade 4. Here we add 
new observations and derive the combined interferomet- 
ric/spectroscopic orbit which agrees with the astrometric 
orbit. 

SIMBAD lists 187 references to date and provides ba- 
sic data on the primary companion, such as spectral type 
G1V, magnitudes V = 5.19, B - V = 0.60, K s = 3.741, 
and near-solar metallicity (see th e full compilation ofstel- 
lar parameter measurements in Soubi ran et al. I [2010). 
The star is located above the Main Sequence in the 
(V, B — V) col or-magnitude diagram (C MP), being obvi- 
ously evolved. iNordstrom et al. " (2004) e stimate t he age 
of 6 Gyr by isochrone fitting. iNielsen fc Close 1 (|20100 
cite indirect age estimates of 4 Gyr and 6.3 Gyr from the 
lithium line strength and rotation, respectively. The star 
was observed by Spitzer and fou nd to have no debris 
disk jTrilling et al. I I2008I ). Yet. iNakaiima fc Morinol 
(|2012f ) claimed recently that k For is young and be- 
longs to the IC 2391 moving group. In calculating the 
kinematic parameters they overlooked the companion 
and used the Hipparcos proper motion (PM), biased by 
the orbit. We calculate the heliocentric Galactic veloc- 
ity components of K For from the center-of-mass PM 
(|Gontcharov et al. l[2000| ) and the new 7-velocity to be 
[U, V, W] = [-19.5, -16.2, -9.6] km s" 1 (U is positive to- 
wards the Galactic center). 

In Section 2 we derive the updated orbit of the outer 
system AB. Then in Section 3 we evaluate the parame- 
ters of the companions and show that the secondary B is 
over-massive, being itself a close binary. This hypothesis 

2 The WDS keeps outdated tradition of assigning discoverer 
codes even when the duplicity was actually discovered before the 
first position measurement, as in the case of k For 
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Figure 1. Combined orbit of k For AB. Left panel: orbit in the plane of the sky, scale in arcseconds. The insert shows the A"- band image 
taken at Gcmini-S on 2012.67. Right panel: the RV curve with early measures plotted as crosses. 

is confirmed in Section 4 by direct spectroscopic detec- 
tion of Ba and Bb. In the last Section 5 we discuss the 
implications of this finding for multiplicity statistics and 
cxo-planet searches. 



2. THE ORBIT OF AB 

So far, spectroscopic, astrometric, and visual orbits 
of k For were determined almost independently of each 
other. We add fresh measures and combine available 
data in the single orbital solution presented in Fig. m 
The orbital elements and their formal errors are obtained 
by the unconstrained least-squares fit with weights in- 
versely proportional to the measurement errors. The 
elements are listed in Table [1] together with the previ- 
ously publi shed orbits. We have no t used the astromet- 
ric data of iGoncharov fc Kivaeva I (|2Q02f) . so the good 
match between spectro-interferometric and astrometric 
orbits speaks to their veracity. This is not the case 
of spectroscopic orbits w hich have a too short period. 
lAbt & Willmarthl (|2006l) mention that they used "some 
older measures" in the orbit derivation but do not spec- 
ify their source. In the future, with better RV and po- 
sitional coverage, the elements will be determined to a 
much higher precision than here, but they will no longer 
be revised beyond the stated errors. The new orbit and 
the H ipparcos parallax of 45.53±0.82mas (|van Leeuwen I 
I2007T ) correspond to the semi- major axis of 11.44 AU and 
to the mass sum of 2.25±0.13M Q (the uncertainty comes 
mostly from the parallax error) . 

Radial ve locities used i n the orbit calculation are listed 
in Table p lEndl et al. I (pOOl did not publish the RVs; 
three points from their Fig. 6 (1993-1998) are neverthe- 
less included in the present solution by assuming that the 
zero velocity in their plot corresponds to 20.36 km s _1 . In 
this case their dat a match well the single measurement in 
1999 publ ished bvlNidever et al. I (|2002f) and the veloci- 
ties from (|Abt fc Willmarth 1 12006). of which we list here 
only the first and the last points. The latest RV mea- 
sures are obtained by the author in 2012 (see Section@|. 
We also use the RVs measures made in Santiago in 1911- 
1916 with the 2-prism spectrometer ([Campbell fc Moorel 



1928). Although the precision of these early data (crosses 
in Fig. [T] right) is low, they help in constraining the or- 
bital period. 

Positional measures of k For AB and their residuals to 
the new orbit are listed in Table [3J The first measure 
is performed with adaptive optics, the remaining data 
come from speckle interferometry at the SOAR telescope 
in Chile. The latest speckle measure was made at SOAR 
on October 29, 2012. The first speckle resolution is given 
a low weight, considering uncertain calibration in that 
early experimental work. 

3. PROPERTIES OF THE COMPONENTS 

The component A (V = 5.19 mag, B — V = 0.60) 
is located ~ 1.5™ above th e Main Sequence in th e 
(M V ,B-V) CMD, as noted bv lNordstrom et al~l (|200l . 
By fitting isochrones, these authors estimated the mass 
of the main star M\ to be between 1.12 and 1.18 M^. 
We ch ecked this against the isochrones of IGirardi et al. I 
( 2000) and found that the mass of A can be as high as 
1.25 Mq if its age is 4 Gyr (it is just leaving the Main Se- 
quence). It is safe to assume that Ma = 1.20 ± 0.05 M©. 

The photometric data are collected in Table |3J The 
magnitude difference in the Stromgren y and Cousins 
/ filters was determined by speckle interferometry at 
SOAR (6 and 5 independent measures respectively); we 
list here the average values and their formal errors and 
assume AV ~ Ay. The derived quantities are listed in 
brackets. 

The flux ratio in the 1.5 4-1.65 /xm band ( which is 
close to H) was estimated bv lLafreniere et al~l (|2007D as 
AH = 2.7 ± 0.2 mag. On 2012 September 2 (2012.6706) 
the binary was resolved with the NICI adaptive optics in- 
strument (|Chun et al. Il2008f ). see Fig. ED As the primary 
companion was heavily saturated, the position measure- 
ment (340.5°, 0'.'342) is not very accurate (it is not in- 
cluded in Table [3]), and no relative photometry could be 
made. However, the flux from the well-resolved compan- 
ion B (after halo subtraction) was estimated by compar- 
ing it to stars HIP 8674 and HIP 12425 observed before 
and after this target at nearly the same airmass. The H 
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Table 1 

Orbits of K For AB 
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Ref. 


(yr) 


(BY) 
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km s 1 


km s _1 




26.5±2 


1983.6±2.3 


0.2±0.1 


0.26±0.01 


133±6 


202±31 


56±4 






a 


21.46±0.09 


1967.26±1.9 


0.30±0.04 






263±12 




3.95±0.26 


18.13±0.66 


b 


26.5 * 


2015.69 


0.20 * 


0.493 


131.6 


275.5 


44.8 






c 


25.81±0.15 


1988.89±0.17 


0.339+0.013 


0.521±0.004 


139.8±1.4 


266.3±1.0 


50.4±0.5 


5.23±0.13 


16.67±0.06 


d 



References: a — astromctric jGoncharov & Kivacva 2002); b 
2012); d - combined, this work 

Table 2 

Radial velocities and residuals of re For A 



JD 


RV 


a 


O-C 


Ref. 


+2400000 


km s — 1 


km s _1 


km s" 1 




19378.67 


18.6 


2.0 


1.08 


a 


19649.90 


19.9 


2.0 


0.58 


a. 


19762.62 


17.4 


2.0 


-2.52 


a 


21152.78 


17.8 


10.0 


-3.69 


a 


49000.0 


21.82 


0.2 


0.04 


b 


50000.0 


20.74 


0.2 


0.00 


b 


50800.0 


19.36 


0.2 


-0.05 


b 


51170.080 


18.824 


0.1 


0.06 


c 


51884.735 


17.71 


0.1 


0.23 


d 
d 


53001.665 


15.46 


0.1 


0.04 


d 


55979.533 


11.53 


0.5 


-0.055 


c 


55983.528 


11.54 


0.5 


-0.052 


c 


56165.830 


12.001 


0.01 


-0.005 


e 


56167.865 


12.019 


0.01 


0.007 


c 


56171.797 


12.032 


0.01 


0.008 


c 


56182.861 


12.048 


0.01 


-0.009 


e 


56194.830 


12.095 


0.01 


0.002 


c 


56200.834 


12.111 


0.01 


-0.001 
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References: 


a - (Campbell & M 


oore 1119281): b 


- lEndl et al. 


2002). zero-point 20.36km e 




— dNidever et 


al. 1 I2002D 


d - i|Ab 


I & Willmarth 


I2006D: e - this 


work. 









Table 3 

Positional measures and residuals of re For AB 
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2005.6348 


267.1 
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0.4 


-2 


a 


2007.8130 
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100 


-1.7 


41 


b 


2008.6282 


294.8 


502 


2 


-0.2 
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2008.6936 


295.3 


501 


2 


-0.3 


-1 


c 


2008.8549 


296.7 


502 


2 


-0.4 





c 


2010.9655 


317.8 


464 


2 


0.6 
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2012.8300 


343.0 


345 
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0.1 
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References: 


a — 


(ILafreniere et al. 1 1200711: 


b - 



2008); 



120101) : d - HHartkopf et al. ll20l3) ; e - this work. 



Table 4 

Photometry of the resolved components 



Parameter 


V 




H 






mag 


mag 


mag 


mag 


t»A+B 


5.19 


4.51 


3.712 


3.741 


m B - m,A 


5.02±0.04 


3.69+0.05 


(2.46) 


(2.14) 


m B 


(10.21) 


(8.23) 


6.3+0.2 


6.0+0.2 



spectroscopic (Abt & Willmarth 2006); c - visual (Hartkopf et al. 

and K s magnitudes of those two stars from 2MASS were 
used to determine the zero point (the actual wavelength 
of the narrow-band filters was 1.587//m and 2.272 //m). 
In this way we obtained crude estimates of the com- 
panion's infrared magnitudes H — 6.0 ± 0.2 mag and 
K s = 6.3 ± 0.2 mag. 

The magnitudes and colors of the companion B are 
thus established reasonably well, (V — I)s = 1.98 ± 0.06 
and (V — K s )b = 4.2 ± 0.2. Standard relations match 
those colors for a dwarf star of ^0.48 M©, but the com- 
ponent B is located about 1.5 m above the Main Sequence 
in the (M v , V - I) and (M V ,V - K) CMDs; its lumi- 
nosity corresponds to a single dwarf of ^0.65 M Q . This 
discrepancy is caused by the binary nature of B. 

Large total mass of the companion B follows from the 
orbital parameters. Motion of the photo-center is charac- 
terized by the semi-major axis a of the astrometric orbit. 
Its ratio to the semi-major axis a of the visual orbit is 
related to the mass ratio q = Mb /Ma and to the light 
ratio r as 



q — r 



a (g+l)(r + l) 



(1) 



In our case we can neglect the companion's light because 
r « 0.01. Th e orbit of AB derived h e re an d the astro- 
metric orbit of iGoncharov k, Kivaeva I (|2002t) correspond 
to 4> = 0.50 and q = (f)/(l — 4>) = 1.0. The companion 
B is therefore as massive as A while being ^100 times 
fainter at optical wavelengths. This finding agrees with 
the large orbital mass sum. 

Another argument for the high companion's mass 
is furnished by radial velocities. Using the formula 



aisini = 0.01375^1^(1 - e 2 ) 



2^.5 



(a in 10 6 km, K x in 
5.53 AU and 



km s , P in days) we obtain a\ 
ai/a = 0.48, leading to q = 0.93. 

The Hipparcos astrometry agrees very well with the 
new orbit if we adopt the companion's masses inferred 
from the photometry, Ma = 1.20 M Q and Mb = 
0.96 M Q , or q = 0.80 and = 0.43. Using this 0, 
we compute the photo-center motion during the Hippar- 
cos mission (mean epoch 1991.25, duration 3.2 yr). Ta- 
ble [5] shows that the orbital proper motion of the photo- 
center A/i and its acceleration fi match the ephemeris 
in both magnitude and direction. In 1991.25 the com- 
panion moved mostly to the South (see the circle in 
Fig. [I] left), the photo-center moved to the North. 
The measured Afi is the difference between the Hippar- 
cos PM of (+197, —5) mas yr -1 and the long-term PM 
of (+ 196, —60) mas yr -1 derived by iGontcharov et al. I 
(2000) from the combination of all ground-based data. 
The Hipparcos astrometry thus confirms the orbit and 
the fact that B is massive. 
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Table 5 

Hipparcos astrometry and the orbit 



Parameter 


Hipparcos 


Orbit 


= 0.43) 




RA Dec 


RA 


Dec 


Afi, mas yr — 1 


1 +55 


-9 


+56 


fi, mas yr -2 


+13.0±2.0 -14.4±1.8 


+ 17.7 


-9.2 



4. SPECTROSCOPY 
4.1. Observations 

Knowing that the companion B must be a close binary, 
we tried to detect this sub-system spectroscopically, de- 
spite its small contribution to the combined light. Opti- 
cal spectra with resolution A/ A A = 80 000 were recorded 
with the fiber-fed CHIRON echelle sp ectrometer in- 
stalle d at the 1.5-m telescope at CTIO (Sc hwab et al. I 
I2012T ) and operated in service mode. The object was 
observed in 2012 six times, on August 25, 27, 31 and 
September 11, 23, 29 (hereafter nights 1 to 6). On each 
visit, two 300-s exposures with the image sheer were 
taken, accompanied by the comparison spectrum of the 
thorium-argon (ThAr) lamp. During this period, the po- 
sition of the ThAr spectrum remained stable on the CCD 
to better than 1 pixel. Moreover, on the first night the 
star HIP 14086 was observed as a RV reference. 

Extracted and wavelength-calibrated spectra were de- 
livered by the pipeline running at the Yale University. 
They contain 59 echelle orders (central wavelengths from 
4605A to 8713A) of 3200 pixels each (1.0768 km s' 1 per 
pixel). As red orders with A > 7000 A suffer from extrac- 
tion and calibration problems (too few thorium lines), 
they are not used here. The maximum intensity in the 
spectra varies from night to night with a total range of 
two times (from 25 to 50 thousand electrons per pixel). 
This corresponds to a S/N from 160 to 220. 

4.2. Radial velocities 



continuum level. The cross-correlation is computed in 
the wavelength range from 4600A to 6500A, thus avoid- 
ing the contamination by telluric lines at longer wave- 
lengths. Regions within ±0.5A of the hydrogen Bahlmcr 
lines are excluded from the mask. The minimum in the 
cross-correlation function (CCF) is approximated by a 
Gaussian curve, its center is taken to be the apparent 
stellar radial velocity, which is then corrected to the 
barycenter of the solar system in the standard way. This 
procedure relies on the wavelength calibration of the re- 
duced spectra. We attempted to refine the velocity zero 
point by cross-correlating red orders with the mask of 
telluric lines, but did not obtain any trustworthy results, 
probably because of the poor wavelength calibration in 
the red. 

Figure [2] shows the RVs derived from 12 individual 
spectra as a function of time (the average RV for each 
night is listed in Table [2]). The linear trend of 2.9 m s" 1 
per day is caused by the orbital motion. The rms scatter 
of RV around this line is 6.8 m s _1 , the rms residual to 
the orbit is 6.7 m s . As measures from the two nightly 
spectra agree well, most of the scatter can be attributed 
to the wavelength calibration based on the ThAr lamp. 
The mean FWHM of the Gaussian fits to the CCF is 
12.96 km s _1 , their equivalent width is 5.59 km s _1 . 

The same procedure applied to HIP 14086 gives an RV 
of +39. 14 km s _1 . According to iNidever et al~l ()2002l ) 
the RV of this star is co nstant at +42.718 km s" 1 , but 
iNordstrom et al. 1 (|2004j ) list a quite different value of 
+38.0km s _1 (while SIMBAD lists the wrong RV of 
—38.4km s _1 ). Apparently, this star is a yet unrecog- 
nized slow spectroscopic binary, unsuitable as a RV stan- 
dard. However, we use our RV measure of HIP 14086 
to derive the RV of k For from the two spectra taken 
in February 2012 during CHIRON tests (those spectra 
lack wavelength calibration). These two measures with 
a somewhat uncertain zero point are given low weight in 
the orbit solution by adopting errors of 0.5km s _1 . 




180 190 
JD - 2456000 

Figure 2. Radial velocity of the component A vs. time. The 
dashed line shows a linear fit with a slope of 2.9 m s — 1 per day. 

Radial velocities of K For were derived by cross- 
correlating the spectra with a binary mask which equals 
one in the lines of the solar spectrum and zero otherwise. 
The mask was constructe d from the digital solar spec- 
trum HHinkle et al. M2Q00rFl binary-clipped at 0.6 of the 

3 NOAO data archive: ftp://ftp.noao.edu/catalogs/arcturusatlas 



4.3. Variability of the hydrogen line profile 

We detect a tiny variability of the Ha profile presum- 
ably caused by the moving emission of Ba and Bb. This 
is achieved by subtracting the template obtained by av- 
eraging all spectra together. Although the spectrometer 
was very stable during these observations, the spectrum 
moved on the CCD by 12 km s _1 owing to the variable 
barycentric correction. The template was constructed by 
shifting the nightly-averaged spectra (with a 3-pixel me- 
dian smoothing of each spectrum to remove the cosmic- 
rays spikes) to match the night 1, normalizing them in in- 
tensity, and adding together. Obviously, the telluric lines 
do not move together with the stellar lines and therefore 
show up in the residuals between the individual spectra 
and the template. The rms residuals after template sub- 
traction are 0.7% in the order 37 containing Ha and 0.5% 
in the blue orders which are free from the telluric lines. 

We see a clear residual signal in the Ha line on 
the nights 1 and 5, while on the remaining nights the 
residuals are close to zero (Fig. [3]). Similar features 
of lower amplitude are also seen in the H/3 line, con- 
firming their reality. The single "emission" peak on 
the night 5 can be approximated by Gaussian curves. 
Their amplitude is 2.1% and 1.3% for Ha and H/? re- 
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Figure 3. Variable features in the H« line on the night 1 (top 
panel) and on the night 5 (bottom panel). In each panel the up- 
per curves show the continuum- normalized spectrum scaled by 0.1, 
with a Gaussian fit to the Ha (dashed line). The lower curves show 
the residuals after template subtraction. 

spectively, and they arc shifted relative to the stellar 
spectrum by +0.28A and +0.23A which corresponds 
to +12. 8km s _1 and +14. 2km s _1 . If the emission 
originated on the slowly-rotating primary companion A 
(Vsini — 4 km s -1 ), it would be centered on the stel- 
lar line rather than shifted. We show below that on 
the night 5 the lines of Ba and Bb were superimposed 
near their center-of-mass velocity which should be dis- 
placed by +9.2 km s _1 relative to the RV of A according 
to the orbit (the masses of A and Ba+Bb being near- 
equal). The double peak on the night 1 also matches 
the expected position of the emission from Ba and Bb. 
The non-detection of emission features on the other four 
nights could possibly be explained by variable activity of 
the M-dwarfs. 

4.4. Retrieval of the Ba and Bb signature by 
cross- correlation 

We detect the signature of absorption lines belonging 
to Ba and Bb by correlating the residual spectra (af- 
ter subtracting the template, see above) with the so- 
lar or Arcturus masks (the latter is analogous to the 
solar binary mask but uses the digital Arcturus spec- 
trum from the same source). Only wavelengths shorter 
than 6500A are used in the correlation. The resulting 
CCF has a weak narrow feature centered on the RV of A 
(Fig. |4]). This feature originates because the resolution of 
the template spectrum is slightly less than the resolution 
of the nightly spectra owing to small residual alignment 
errors in the template creation. The sign and intensity of 
this feature depend on the relative degree of smoothing 



TRIPLE RADIO-STAR 5 

applied to the spectra and the template. 
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Figure 4. Correlation functions of residuals with the Arcturus 
mask for 6 nights (from bottom up, displaced by 0.0015 per night). 
The horizontal axis corresponds to the heliocentric RV. The upper 
curve is correlation with the template to which a fake companion 
of 0.5% intensity and +100 km s — 1 RV shift is added. The dashed 
and dotted lines show suggested positions of the dips corresponding 
to Ba and Bb components respectively moving on a 3.7-d orbit. 
Numbers on the left are Julian dates, numbers on the right are the 
nights. 

Apart from the central feature, we see additional dips 
with variable position and intensity, presumably pro- 
duced by Ba and Bb. Gaussian approximations of the 
strongest dips have FWHM between 15 and 23 km s _1 
and equivalent width of about 0.0017 km s _1 . The star 
gives a CCF with equivalent width of 3.54 km s _1 for 
the Arcturus mask, so the dip ratio corresponds to the 
flux difference on the order of 5.8 m (we neglect here the 
different match of the Arcturus mask with the late-type 
spectrum of Ba on one hand and with the GOV spectrum 
of A on the other hand). A fake companion of 0.5% rela- 
tive intensity (5.75 m fainter than the primary) is readily 
detectable by this method (see the upper curve in Fig.@|. 

The curves in Fig. |4] give an impression of secondary 
lines moving with an amplitude of ^80 km s . We 
found a tentative circular orbit of Ba,Bb with a period 
of 3.666 d, semi- amplitude K\ = K% = 83 km s , 7- 
velocity +26 km s , and initial epoch (RV maximum) 
JD 2456166.45. The data at hand are not sufficient for 
the orbit determination and the above elements are a 
guess only. The dashed and dotted lines in Fig. |4] are 
Gaussian dips corresponding to the Ba and Bb compo- 
nents respectively and positioned according to this ten- 
tative orbit (not fitted to the CCF). The assumed am- 
plitudes of both Gaussians are -0.0007, their FWHM 
is 16 km s _1 . On the night 1 we clearly see both dips. 
On the nights 3 and 4 the CCFs of residuals look almost 
identical, the dips are near the maximum separation. Re- 
member that the template includes the average spectrum 
of Ba and Bb, so their lines near maximum separation 
(the most frequent situation) are partially subtracted, 
leaving only small residuals. On the night 5 the dips 
overlap near the 7-velocity. 

The reality of the detection of Ba and Bb was checked 
by correlating the residuals with the synthetic spec- 
trum of a late- type dwarf. We used the R = 500 000 
spectrum from (fBertone et al. I l2b08fl with parameters 
T e = 4000 K, log 5 = 5.0, and solar metallicity. The re- 

4 http : //www. inaoep.mx/- Omodelos/bluered/documentat ion. html 
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Figure 5. Correlation functions of residuals with the M-star syn- 
thetic spectrum for 6 nights (from bottom up), similar to Fig. [4] 

suiting CCF should have a peak corresponding to the 
velocity of the secondary. Indeed, we see in Fig. [5] posi- 
tive features at approximately same velocities as the dips 
in Fig. [4] To highlight these features, we again over-plot 
fiducial Gaussians with an amplitude of +0.00007 and 
velocities predicted by the tentative Ba,Bb orbit. 

The case of the component B being a binary is very 
strong. Its high mass, position in the CMD, variable 
emission, and moving absorption lines all confirm this 
hypothesis. The RV semi-amplitude in a spectroscopic 
binary with a circular edge-on orbit of period P days is 
A a = 213P" :L / 3 M 2 (Mi + M 2 y 2/3 km s _1 , or 69km s" 1 
for P = 3.7 d, assuming M x = M 2 = 0.5M Q . The RV 
amplitude of our tentative Ba,Bb orbit is therefore in 
qualitative agreement with its period. 

5. DISCUSSION 

The triple system k For consists of the main slightly 
evolved (age 4-6 Gyr) component A of 1.2 Mq, and a 
pair of M-dwarfs of ~0.5M Q each on a tight orbit with 
a tentative period of 3.7 d. The masses of A and Ba+Bb 
are nearly equal. 

Further spectroscopic observations of k For are needed 
to confirm the tentative period of Ba,Bb and to deter- 
mine its orbit and mass. High-resolution spectroscopy in 
the near-infrared where the contrast of Ba,Bb is more fa- 
vorab le will help in this endeavor (e.g. iBender fc SimorTI 
2008). Given the short period, eclipses in Ba,Bb are 
likely, but their photometric detection presents a chal- 
lenge because of the small (sub-percent) amplitude in 
the combined visible light. Precise photometry will even- 
tually detect flares on the active M dwarfs. Long-term 
monitoring of both orbits will lead to the accurate mea- 
surement of the distance and the masses of all three stars. 

Until now, the triple nature of k For was not recog- 
nized despite its closen ess to the Sun. The survey of 
iRaghavan et al. 1 (2010) lists it only as a binary. The 
question is how many more such systems are we miss- 
ing? This relates to the ways of discovering sub-systems 
in the faint secondary companions. 

The binarity of B was revealed by its unusually high 
mass. Over-massive companions can be detected by as- 
trometry quite easily, comparing the photo-center axis a 
with the estimated full axis a (see eq. [T]), if the astro- 
metric orbit and parallax are known. However, reliable 
astrometric orbits are rare. Acceleration measured by 
Hipparcos can also be used to estimate the mass ratio 
in binaries with reliably known visual or astrometric or- 



bits. On the other hand, an excess of the total mass 
in orbital visual binaries with known parallax is a less 
promising way of finding such sub-systems because of 
the typically large uncertainties in both measured and 
modeled masses. Over-massive (binary) secondaries also 
can be detected by a high minimum mass inferred from 
a spectroscopic orb it, as e.g. in the case of the quadruple 
syste m HD 27638 (|Tokovinin fc Gorvnva 1120011 : iTorres I 
2006). Over -massive but invisible companions could be 
white dwarfs; distinguishing them from tight pairs of red 
dwarfs usually requires some sort of photometry. When 
the resolved photometry of the secondary companion in 
two colors is available (e.g. from speckle interferometry) , 
it can be placed on the CMD and detected as a binary by 
an excess in luminosity. Last but not least, spatial reso- 
lution of the sub-systems is the most direct and powerful 
method; if the period of Ba,Bb were longer than ~0.5yr 
(semi-major axis > 30 mas), it would have been resolved 
at SOAR. 

Statistical modeling done by iTokovinin et al. 1 ()2012f ) 
hints that 10% to 20% of nearby solar-type astrometric 
binaries found by Hipparcos could have "massive" secon- 
daries (sub-systems or white dwarfs). This conclusion, 
however, depends on several assumptions in the model. 
The detection techniques outlined above can be applied 
systematically to nearby binaries to better constrain the 
fraction of massive secondaries without making any as- 
sumptions. If sub-systems in the secondary companions 
are indeed as frequent as in the primaries, the known 
number of hierarchical triples will be nearly doubled, 
changing the multiplicity statistics dramatically. 

Origins of multiple stars are still actively researched 
and debated, and different theories are tested by their 
predictions about multiplicity. For example, chaotic N- 
body dynamics rarely produces pairs of low-mass stars 
revolving around more massive primaries; rather, the 
primary itself is likely to end up in a close pair while 
the least massive body is ejected into a distant orbit 
around it (Dclgado -Donate et al. Il2003l ). Fragmentation 
of rotating cores makes the opposite prediction where 
the sub-systems in primary or secondary components are 
equally likely to form, being a natural sink of the angu- 
lar momentum. In this paradigm, stars are often born 
as quadruples, but some inner sub-systems subsequently 
merge. Originally k For could have consisted of two pairs 
of similar 0.5-Mq stars. One pair had merged long time 
ago and became the present-day primary, while the other 
pair Ba,Bb is still here. 

Binary secondary companions such as the one of k For 
have implications for the search of exo-planets. If the 
period of Ba,Bb were few months or years, the weak 
secondary lines would always blend with the lines of A 
and create a small periodic RV signal that could be mis- 
taken for an exo- planet signature, as e.g. in HD 19994 
(|Roll et al. Il2010h . 
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